Stochastic Quantization of Scalar Fields 
in Einstein and Rindler Spacetime 



G. Menezes ^ and N. F. Svaiter ^ 

Centro Brasileiro de Pesquisas Fi'sicas-CBPF, 
Rua Dr. Xavier Sigaud 150, 
Rio de Janeiro, RJ, 22290-180, Brazil 



Abstract 

We consider the stochastic quantization method for scalar fields defined in a curved manifold 
and also in a fiat space-time with event horizon. The two-point function associated to a massive 
self-interacting scalar field is evaluated, up to the first order level in the coupling constant A, 
for the case of an Einstein and also a Rindler Euclidean metric, respectively. Its value for the 
asymptotic limit of the Markov parameter r — oo is exhibited. The divergences therein are taken 
care of by employing a covariant stochastic regularization, where all the symmetries of the original 
theory are preserved. 
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1 Introduction 



In the last century Parisi and Wu introduced the stochastic quantization method [1]. The 
main idea of the stochastic quantization is that a rf-dimensional quantum system is equivalent 
to a ((i + l)-dimensional classical system including fluctuations effects into considerations. Some 
of the most important papers in the subject can be found in Ref. [2]. A brief introduction to 
the stochastic quantization can be found in the Refs. [3] and [4] and a complete review of the 
stochastic quantization is given in the Ref. [5]. 

This program of stochastic quantization and the stochastic regularization was carried out for 
generic flelds deflned in flat, Euclidean manifolds. It is important to observe that, since the 
stochastic regularization is not an action regularization, this method preserves all the symmetries 
of the theory under study. In the development of this program some authors applied this method 
to linearized Euclidean gravity and also non-linearized gravity [6] [7] [8] [9] [10]. It is worth 
pointing out here that we cannot expect that the method will improve the pcrburbativc non- 
renormalizability of quantum gravity. Working along such lines, we may observe that the study 
of a situation which lies between these two extremes is missing. Therefore, before the attempt to 
implement the program in non-renormalizable models, a consistent logical step is to discuss an 
intermediate situation between fields in fiat spacetime and quantum gravity, i.e., the semiclassical 
theory, which corresponds to the one-loop approximation to the totally quantized theory. It is 
tempting to think that many unsolved problems, as the loss of information in black hole formation 
and the entropy of black holes, can gain new insights by introducing the stochastic quantization. 

Although there are some papers discussing the stochastic quantization in Minkowski spacetime 
[11] [12], the stochastic quantization was originally introduced as an alternative quantization 
method for fields defined in Euclidean space. Therefore our aim in this article is to discuss the 
stochastic quantization of scalar fields defined in a static curved manifold without event horizon 
and also in a flat manifold with event horizon i.e., to analytically continue metrics that generate 
real actions. We would like to stress that the stochastic quantization is quite different from the 
other quantization methods, therefore it can reveal new structural elements of a theory which so 
far have gone unnoticed. 

The method of stochastic quantization in flat spacetime with trivial topology can be summa- 
rized by the following steps. First, starting from a field defined in Minkowski spacetime, after 
analytic continuation to imaginary time, the Euclidean counterpart, i.e., the field defined in an 
Euclidean space, is obtained. Second, it is introduced a monotonically crescent Markov parameter, 
called in the literature "fictitious time" and also a random noise field ?7(r, x), which simulates the 
coupling between the classical system and a heat reservoir. It is assumed that the fields defined at 
the beginning in a (i-dimensional Euclidean space also depends on the Markov parameter, there- 
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fore the field and a random noise field are defined in a (rf + l)-dimcnsional manifold. One starts 
with the system out of equilibrium at an arbitrary initial state. It is then forced into equilibrium 
assuming that its evolution is governed by a Markovian Langevin equation with a white random 
noise field [13] [14] [15]. In fact, this evolution is described by a process which is stationary, 
Gaussian and Markovian. Finally, the n-point correlation functions of the theory in the {d + 1)- 
dimensional space are defined by performing averages over the random noise field with a Gaussian 
distribution, that is, performing tlic stochastic averages {'^{ri,xi)ip{T2,X2)---'^{Tn,Xn))ri- The n- 
point Schwinger functions of the Euclidean d-dimensional theory are obtained evaluating these 
n-point stochastic averages ( v?(ri, Xi)v?(r2, a;2)...</?(r„, x„) when the Markov parameter goes to 
infinity (r oo), and the equilibrium is reached. This can be proved in different ways for the 
particular case of Euclidean scalar field theory. One can use, for instance, the Fokker-Planck 
equation [16] [17] associated with the equations describing the stochastic dynamic of the system. 
A diagrammatical technique [18] has also been used to prove such equivalence. 

The original method proposed by Parisi and Wu was extended to include theories with fermions 
[19] [20] [21]. The first question that appears in this context is if make sense the Brownian problem 
with anticommutating numbers. It can be shown that, for massless fermionic fields, there will 
not be a convergence factor after integrating the Markovian Langevin equation. Therefore the 
equilibrium is not reached. One way of avoiding this problem is to introduce a kernel in the 
Langevin equation describing the evolution of two Grassmannian fields. 

Usually, the Parisi- Wu scheme of quantization applied to bosonic and fermionic fields converges 
towards non-regularized theories. In order to obtain regularized stochastic diagrams, the original 
stochastic process proposed by Parisi and Wu must be modified. One way is to replace the 
Markovian process by a non-Markovian one [22] [23]. This can be done introducing a colored noise 
in the Einstein relations [24] [25] [26] [27]. Recently, the stochastic quantization of a topological 
Chern-Simons theory and the self-interacting X(f^ theory was investigated by the authors, using 
the Parisi- Wu scheme with a non-Markovian Langevin equation [28] [29]. 

Motivated by the problems that appear when one tries to implement this regularization scheme 
in non-abelian gauge theories with Zwanziger's gauge fixing [30] [31] [32], Bern et al [33] derived 
a truly non-perturbative regularization procedure in QCD, still in the Markovian scenario. They 
introduced a modified Langevin equation with a regulator multiplying the noise. Therefore, the 
original Einstein relations and consequently the noise probability distribution are maintained. 
This continuum regularization program was very successful in the invariant non-perturbative reg- 
ularization of all quantum field theories [34] [35] [36], including gauge theory, gauge theory with 
fermions and gravity. In few words, the idea of any stochastic regularization is to start from an 
interacting theory, then construct Langevin tree-graphs, where each leg of which ends in a regu- 
larized noise factor. Since it is possible to obtain the loops of the theory by contracting the noise 
factors, one ends up with a theory where every closed loop contains at least some power of the 
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regulator. With this modification, one can show that the system converges towards a regularized 
theory. The next step to construct a finite theory is to use for example a minimal-subtraction 
scheme in which the ultraviolet divergent contributions are eliminated. 

The aim of this paper is to implement the stochastic quantization for the self-interacting A<^^ 
theory in a Riemannian manifold and also in a flat manifold with event horizon, being more 
speciflc, the Einstein [37] and the Rindler spacetime [38] [39] [40] [41], where positive-deflnite 
forms can be obtained which generate real actions in the Parisi-Wu Langevin equation. The only 
article that has some overlap with our work is the Huang paper [42] . This author have used the 
stochastic quantization to find the regularized stress-energy tensor associated to a scalar field in an 
inhomogeneous spacetime. We would like to stress that although Brownian motion on a manifold 
has been studied by the mathematicians, there are a few papers studying Brownian motion in 
general manifolds. See for example the Ref. [43]. 

This quantization method differs from the others, the canonical and the path integral field 
quantization, based in the Hamiltonian and the Lagrangian, respectively, in many aspects. The 
method starts from a classical equation of motion, but not from Hamiltonian or Lagrangian, and 
consequently can be used to quantize dynamical systems without canonical formalism. Further- 
more, it is useful in situations where the others methods lead to difficult problems and can bring 
us new important results. 

For example, a quite important point in a regularization procedure is that it must preserve all 
the symmetries of the unregularized Lagrangian. Many authors have been stressed that a priori 
we can not expect that a regularization independent proof of the renormalization of theories in a 
curved background exists. The presence of the Markov parameter as an extra dimension lead us 
to a regularization scheme, which preserves all the symmetries of the theory under study. Since 
the stochastic regularization is not an action regularization, it may be a way to construct such 
proof. As a starting point of this program, we should calculate the two-point function up to the 
first order level in the coupling constant A and apply the continuum stochastic regularization. Our 
results are to be compared with the usual ones in the literature. 

The organization of the paper is the following: in section H we discuss the stochastic quantiza- 
tion for the {\(f^)d scalar theory in a d-dimensional Euclidean manifold. In section HI we use the 
stochastic quantization and the stochastic regularization to obtain the two-point Schwinger func- 
tion in the one-loop approximation in the Einstein manifold using the stochastic quantization. In 
section IV we repeat the method for the case of the Rindler manifold, a non-simple connected man- 
ifold in the Euclidean version. Conclusions are given in the section V. In the appendix we sketch 
the general formalism to construct the maximal analytic extension of the vacuum Schwarzschild 
solution. In this paper we use h — c — ks — G— 1. 
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2 Stochastic quantization for the {X(p^)d scalar theory: the 
Euclidean case 

In this section, wc give a brief survey for the case of self-interacting scalar fields, implementing 
the stochastic quantization and the continuum stochastic regularization theory up to the one-loop 
level. Let us consider a neutral scalar field with a {Xcp'^) self-interaction. The Euclidean action 
that usually describes a free scalar field is 

SoM = I d'x (^^{d^r + Imlcp'ix)^ , (1) 

and the interacting part, defined by the non-Gaussian contribution, is 

SiM = jd'xjyix). (2) 

The simplest starting point of the stochastic quantization to obtain the Euclidean field theory 
is a Markovian Langevin equation. Assume a flat Euclidean (/-dimensional manifold, where we are 
choosing periodic boundary conditions for a scalar field and also a random noise. In other words, 
they are defined in a d-torus Vl = T'^. To implement the stochastic quantization we supplement the 
scalar field ip{x) and the random noise 'r]{x) with an extra coordinate r, the Markov parameter, 
such that (p{x) — > <f{T,x) and 77(0;) — > r]{T,x). Therefore, the fields and the random noise are 
defined in a domain: T ^ x R^~^\ Let us consider that this dynamical system is out of equilibrium, 
being described by the following equation of evolution: 

^V?(r,x) = —^^Ux)=^(t,x) + Vir^x), (3) 

where r is a Markov parameter, 77 (r, x) is a random noise field and Sq is the usual free Euclidean 
action defined in Eq.(l). For a free scalar field, the Langevin equation reads 

d 

— (/:>(r, x) = -(-A + ml )ip{T, x) + r]{T, x), (4) 

where A is the (i-dimensional Laplace operator. The Eq.(4) describes a Ornstein-Uhlenbeck process 
and we are assuming the Einstein relations, that is: 

(r/(r,x)), = 0, (5) 

and for the two-point correlation function associated with the random noise field 

{v{t,x)v{t',x')), =26{t-t')6\x-x'), (6) 
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where ( ...)ri means stochastic averages. The above equation defines a delta-correlated random 
process. In a generic way, the stochastic average for any functional of ip given by F[(p] is defined 

by 



J[dr]]F[(p]ex.p 




x) 


J [dr]] exp 




iJd'^xJdTTj'^iT, x) 





(7) 



Let us define the retarded Green function for the diffusion problem that we call G{t — t',x — x'). 
The retarded Green function satisfies G{t — t' ,x — x') — Q \l t — t' < Q and also 



d_ 
dr 



G{t - r', X - x') = 5\x - x')5{t - t'). 



(8) 



Using the retarded Green function and the initial condition (^(r, x)\t=q — 0, the solution for Eq.(4) 

reads ^ 

(p{T,x)= [ dr' [ d'^x'G{T-T',x-x')ri{T',x'). (9) 
JO Jn 

Let us define the Fourier transforms for the field and the noise given by (/?(r, k) and r7(r, k). We 
have respectively 

" d'^xe-'^''ip{T,x), (10) 



and 



r/(T, k) 



J d^xe-'^''r]{T,x). 



(11) 



Substituting Eq.(lO) in Eq.(l), the free action for the scalar field in the {d+ l)-dimensional space 
writing in terms of the Fourier coefficients reads 



So[ip{k)] \^p{k)=Lp{T,k) 



d'^kip{T, k){k'^ + ml)ip{T, k). 



(12) 



Substituting Eq.(lO) and Eq.(ll) in Eq.(4) we have that each Fourier coefficient satisfies a Langevin 
equation given by 

|-(^(t, k) ^-{k^ + mlMr, k) + 7y(T, k). (13) 

In the Langevin equation the particle is subject to a fluctuating force (representing a stochastic 
environment), where its average properties are presumed to be known and also the friction force. 
Note that the "friction coefficient" in the Eq.(13) is given by (/c^ + mp). 
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The solution for Eq.(13) reads 

(/p(t, k) = exp [-(e + mDr) (^(0, k) + £ dr' exp (-{k^ + m^)(r - r')) k). (14) 

Using the Einstein relation, we get that the Fourier coefficients for the random noise satisfies 

(r7(T,fc)), = (15) 

and 

( 77(t, ^)77(t', k') ), = 2S{t - r')5\k + k'). (16) 
It is possible to show that ( (^(r, k)(p{T', k') )r,|r=r' = D{k, k'; r, r') is given by: 

D{k; T, r) = (27r)'^5'^(^ + /c')^^^^^ (l - exp (-2t(A;2 + m^))) . (17) 

where we assume r = r'. 

Now let us analyze the stochastic quantization for the {Xip'^)d self-interaction scalar theory. In 
this case the Langevin equation reads 

d A 

— <^(r, x) = -(-A + ml )if{T, x) - ^<^^(t, x) + r]{T, x). (18) 

The two-point correlation function associated with the random field is given by the Einstein 
relations, while the other connected correlation functions vanish, i.e., 

{r]{Ti, Xi)r]{T2, X2) .■.'n{r2k-l, X2k-l) )r, = 0, (19) 

and also 

{ri{Ti,Xi)...r]{T2k,X2k) )ri = Y. ivin, Xi)ri{T2, X2) )n{viniXk)'n{Ti,Xi) )^..., (20) 

where the sum is to be taken over all the different ways in which the 2k labels can be divided 
into k parts, i.e., into k pairs. Performing Gaussian averages over the white random noise, it is 
possible to prove the important formulae 

T / / ^ r ^ f N\ J[d(p]^{xi)^{x2)...^{Xn) 6'^^'^^ 

^lim ( (p{Ti, xi)(p{T2, X2)...ip{rn, Xn) )^ = /[(^(^] e-^(^) ' ^ ^ 

where S{ip) = SQ{ip) + Si{ip) is the d-dimensional action. This result leads us to consider the 
Euchdean path integral measure a stationary distribution of a stochastic process. Note that the 
solution of the Langevin equation needs a given initial condition. As for example 

(p{T,x)\r=o = V^o{x). (22) 
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Let us use the Langevin equation to perturbatively solve the interacting field theory. One way 
to handle the Eq.(18) is with the method of Green's functions. We defined the retarded Green 
function for the diffusion problem in the Eq.(8). Let us assume that the coupling constant is a 
small quantity. Therefore to solve the Langevin equation in the case of a interacting theory we 
use a perturbative series in A. Therefore we can write 

(fir, x) = (r, x) + Xifi^^^ (r, x) + AV^'^ (t, x) + ... (23) 

Substituting the Eq.(23) in the Eq.(18), and if we equate terms of equal power in A, the resulting 
equations are 

ip^'\T,x)^ri{T,x), (24) 

^W(T,x) = -i(^W(r,x))', (25) 

and so on. Using the retarded Green function and assuming that (p^'^^r, x)\t-=o = 0, W q, the 
solution to the first equation given by Eq. (24) can be written formally as 

^(0) (r, x) = r dr' I GiT-T\x- x'Mt' , x') . (26) 

Jo JQ 

The second equation given by Eq.(25) can also be solved using the above result. We obtain 

(^^^^(r, x) = —7^ / dri [ d'^xiGlr — Ti,x — xi) 
3! Jo Jn 

(^1^' dr' l^d'^x'G{Ti-T',xi-x')'n{T',x')^ . (27) 

We have seen that we can generate all the tree diagrams with the noise field contributions. We 
can also consider the ?7.-point correlation function ( (/?(ri, a;i)(/?(r2, X2)...<^(r„, a;„) ),,. Substituting 
the above results in the n-point correlation function, and taking the random averages over the white 
noise field using the Wick-decomposition property defined by Eq.(20) we generate the stochastic 
diagrams. Each of these stochastic diagrams has the form of a Fcynman diagram, apart from the 
fact that we have to take into account that we are joining together two white random noise fields 
many times. Besides, the rules to obtain the algebraic values of the stochastic diagrams are similar 
to the usual Feynman rules. 

As simple examples let us show how to derive the two-point function in the zeroth order 
( (/^(ri, Xi)(/9(r2, X2) and also the first order correction to the scalar two-point-function given by 
{(P{ti, xi)ip{T2, X2) )!^^. Using the Eq.(9) and the Einstein relations we have 

/•min(Ti,T2) r 

{ ip{n,x,)ip{T2, X2) )(°' = 2 dr' d^x' Gin - r', x, - x') G{t2 - r' , X2 - x'). (28) 



d_ 



+ (- 



-A^ + ml 



d_ 



+ (-A^ + ml ) 
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For the first order correction we get: 

= -A( J dX^ j dX^{G{X^ - X^)G{X2 - X3) + - X^)G{X2 - X^)^ 
v{Xs)(^J dX,G{X,-X,)r]{X,)^ (29) 

where, for simphcity, we have introduced a compact notation: 

£dTj^d'^x = J dX, (30) 

and also (/?(r, x) = ifi{X) and finally ^^(r, x) = 'r]{X). 

The process can be repeated and therefore the stochastic quantization can be used as an 
alternative approach to describe scalar quantum fields. Therefore, the two-point function up to 
the first order level in the coupling constant A is given by 

(^(Ti,Xi)^(r2,X2))« = (a) + (6) + (c), (31) 
where (a) is the zero order two-point function and (b) and (c) are given, respectively, by: 

W = -^5'(^i + ^2) / d^'k dTG{h;n-T)D{k;T,T)D{k2;T2,T), (32) 

(c)^-^5\ki + k2) j d'^k£' dr G{k2;T2 - T)D{k;T,r)D{ki;n,T). (33) 

These are the contributions in first order. A simple computation shows that we recover the correct 
equilibrium result at equal asymptotic Markov parameters (ti = r2 — > 00): 

(^')ln=..^oo = S'{k, + ^2) + ^2) ^kl + kl+2ml) I "^'^k^lmlY ^^^^ 

Obtaining the Schwinger functions in the asymptotic limit does not guarantee that we gain a 
finite physical theory. The next step is to implement a suitable regularization scheme. A crucial 
point to find a satisfactory regularization scheme is to use one that preserves the symmetries of 
the original model. The presence of the Markov parameter as an extra dimension lead us to a new 
regularization scheme, the stochastic regularization method, which preserves all the symmetries 
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of the theory under study. Therefore, let us implement a continuum regularization procedure [36]. 
We begin with a regularized Markovian Parisi-Wu Langevin system: 

d 6 S f 

—ip{T,x) = -j^^\^(x)=^{t,x) + J d'^yRxy{^)r]{T,y). (35) 

The Einstein relations, given by Eqs. (5) and (6), are maintained. The regulator -R(A) that 
multiplies the noise is a function of the Laplacian: 

A,, = j d''z{d^U{d^),y, (36) 

where 

{d,)xy^&;,5\x-y). (37) 
We will be working with a heat kernel regulator with the form: 

i?(A;A) = exp(A), (38) 

where A is a parameter introduced to regularize the theory. The basic restrictions on the form of 
this heat kernel regulator are: 

i?(A;A)|A_oo = l, (39) 

or 

i?,,(A;A)|A_,^ = 5'^(x-y), (40) 

which guarantees that the regularized process given by Eq.(35) reduces to the formal process given 
by Eq.(3) in the formal regulator limit A oo. 

With this modification in the Langevin equation, it is possible to show that all the contributions 
to the n-point function at all orders in the coupling constant A are finite. For instance, the 
contribution to the two-point function at the one- loop level given by Eq.(32) is rewritten as: 

(^>)ln=r.->oo = -^^"(^1 + ^2) ,,,^^^ ,J , / d'^kj-^^^, (41) 

V /in r2 oo 2 ^ ' ^2^^ 2ml) J {k^ + mlY ^ ' 

where is the Fourier transform of the regulator, i.e., 

i?,(A) = i?(A;A)U=-fe2. (42) 

Now we use this method to discuss the quantization of scalar theories with self-interaction in 
a curved spacetime without event horizon and in a flat manifold with event horizon. Being more 
specific, we are interested to investigate the A(/9'^ theory in the Einstein and Rindler spacetime, 
respectively. 
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3 Stochastic quantization for the {X(p^)d scalar theory: the 



Einstein case 

The aim of this section is to implement the stochastic quantization and the stochastic regu- 
larization for the self-interacting \(^'^ theory in the one-loop level in the Einstein spacetime. Let 
us consider a manifold that admit a non- vanishing timelike Killing vector field X. If one can 
always introduce coordinates t = x^, x^, locally such that X = and the components of 
the metric tensor are independent of t, is stationary. If further the distribution X-^ of 3-planes 
orthogonal to X is intcgrable, then is static. Each integral curve of the Killing vector vector 
field X = ^ is a world line of an possible observer. Since X = generates isometrics, the 
3-planes X are invariant under these isometrics. For static manifold, it is possible to perform 
a Wick rotation, i.e., analytically extend the pseudo-Riemannian manifold to the Riemannian 
domain without problem. Therefore for static spacetime the implementation of the stochastic 
quantization is straightforward. 

In the previous section, we have been working in an Euclidean space R'^ X R^+\ where 
is the usual Euclidean space and R^'^^ is the Markov sector. Now let us generahze this to a 
more complicated case, i.e., let us work in a general (and, for the time being, static) Riemannian 
manifold M. In other words, we will consider a classical field theory defined in a M x 
manifold coupled with a heat reservoir. Therefore, in a four-dimensional curved manifold, the 
Parisi-Wu Langevin equation for scalar fields reads 

-^(r, x) = -— ^|,(.)=,(.,.) + v(r, x), (43) 

where g — det gf^u and the classical Euclidean action 5*0 is given by 

So^Jd'x^C. (44) 

In the above equation C is given by 

J^^^g^^u d^<pd,ip + ^ (m^ + ^R)^\ (45) 

Note that we introduce a coupling between the scalar field and the gravitational field represented 
by the term ^Ri^^, where C a, numerical factor and R is the Ricci scalar curvature. The random 
noise field r]{T,x) obeys the following generahzed Einstein relations: 

{r}{T,x)), = 0, (46) 



10 



and 

{rj{T,x)ri{T',x'))r, = -jL=S\x-x')diT-T'). (47) 
Substituting Eq.(44) and Eq.(45) in the Langevin equation given by Eq.(43), we get: 

^^(t,x) = -(^-A + m^ + ^B}jip{T,x)+rj{r,x), (48) 
where A is the four-dimensional Laplace-Beltrami operator defined by: 

= (49) 

V denoting the covariant derivative. 

To proceed, as in the flat situation, let us introduce the retarded Green function for the diffusion 
problem G{t — r', x, x'), which obeys: 



d_ 



G{t - t', X, x') = — 5\x - x')5{t - t'), (50) 
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if r - r' > 0, and G{t - r', x, x') = ff r - r' < 0. 

Using the retarded Green function and the initial condition (/?(r, x)\r=o = 0, a formal solution 
to Eq.(48) reads 

(^(r, x) = r dr' [ d^x' \f^)GiT -t',x- x'Wt' , x'). (51) 
Jo Jn 

We see that, in order to solve Eq.(48) we need to invert the differential operator given by Eq.(49). 
Depending on the metric g^i, defined on the manifold, this could be a difficult task. In the 
Euclidean case this is easily solved, since we decompose the fields in Fourier modes (plane waves), 
which are solution to the Klein-Gordon equation in a flat metric; each Fourier mode, as we have 
seen, obeys a Langevin equation. So, if the Klein-Gordon equation in a general curved manifold 
allows us to obtain solutions which can be decomposed in modes separated in a similar manner 
to the flat case, we may as well decompose our flelds in these modes in a way that each of them, 
again, will obey a Langevin equation. Then, Eq.(48) can be solved in a simple and direct way. 
Of course, this procedure will help us if we admit manifolds with global Killing vectors, which, 
unfortunately, is not always the case. So admitting that this is the case, we may write the mode 
decompositions as: 

(f{T,x)= / djl{k)(fk{r)uk{x), (52) 
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and 

r]{T,x)= I djl{k)r]k{T)uk{x), (53) 



where the measure ^{k) depends on the metric we are interested in. For instance, in the flat case, 
we have that in a four dimensional space = d'^k and the modes Uk{x) are given by: 

u,(x) = ^e^^^ (54) 

So we see that, in the flat case, Eqs. (52) and (53) reduce to the Fourier decomposition given by 
Eqs.(lO) and (11), respectively. 

Now let us apply this formalism to the Euclidean Einstein manifold. The general Robertson- 
Walker line element is given by [44] [45]: 

3 

ds'^ = -dt"^ + a^{t) ^ hijdx'dx^, (55) 

where: 

3 _^ 

^hijdx'dx^ = (l-Kr^) dr^ + r^de^ + sin^ (j)d(l)^) 

^dx' + fix)ide' + Sin' cl)d<^'), (56) 



and, for the Einstein manifold we have {K = 1), /(%) = r = sinx, < % < 2n. 

The Eq.(56) gives the line element on the spatial sections in the pseudo-Riemannian case 
which are hyperbolic, flat or closed depending on whether K — —1,0,1, respectively. Writing 
C{p) — a^^t), with the conformal time parameter p given by 



P = 

the line element can be recast in the form: 



y* dt'a-\t'), (57) 



ds' - C(p) (-dp' + ^ hijdx'dxA . (58) 



In the simplest case, namely static space-times with C{p) = C = a' = constant, the scalar 
curvature is given by i? = After the Wick rotation, we are working with the Euchdean 
Einstein manifold. In this case, we may decompose the modes as: 

Uk{x) = C-'/'Xr{x)e''^P, (59) 
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with X = (r, 6, 0) or 6, 0) and Xg being a solution of 

A(^)X,- = -(g^ - K) X,-. (60) 
In Eq.(60), A*^^) is the Laplacian associated with the spatial metric hij: 

A = h-^l^di{h}l^h''dj). (61) 
Further, the functions Xj^ are normalized such that: 

j d^x X^{x)X*g{x) = 5{k, k'), (62) 

where S{k, k') is the 5-function with respect to the measure ji: 

j dfi{k')f{k)5{k,k') = f{k'). (63) 
The eigenfunctions X^^ of the three-dimensional Laplacian are, for K — 1: 

x,-(f) = n+(x)W^», (64) 

with k = {q,J,M), M = -~J,-J + 1,...,J, J = 0,l,...,q- 1 and g = 1,2,.... The Ymj are the 
usual spherical harmonics. The functions II^j can be defined from 



-1/2 / ^ xJ+1 

^^^^^^ Wcoshx i (^^) 



by replacing q by —iq and x by —ix in the latter. 

With these definitions, the measure jl{k) for the Einstein universe is defined as follows: 



dm = ^ I dkpY. ■ (66) 

2^ q,J,M 



So, inserting the mode decomposition given by Eq.(59) in Eq.(48) we have that each mode coeffi- 
cient satisfy the Langevin equation given by 

-^V^kir) = ~{k'' + n^)<pk{r) + %(r), (67) 

where k"^ = q^ + k^ and /i^ = C'm'^ + {6^ — l)K and K — 1. For simplicity, we redefine ^{k'^ + ii^) 
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The solution for Eq.(67), with the initial condition (pk{T)\T=o = 0, reads: 

ipkir)^ rdT'Gk{T,T')rj{T',k), (68) 
Jo 

where 

Gfe(r, r') = exp (-(P + fi'){T - r')) ^(r - r') (69) 

is the retarded Green function for the diffusion problem. 

Using the Einstein relations, we get that the mode coefficients for the random noise satisfies 

(%(r)), = (70) 

and 

{Vk{r)Vk'{r') ), = 25(t - T')S\k, k'), (71) 
where 5^(A;, k') = 5{kp + k'pjSqq'Sjj'SuM'- 

The two-point function Dk{T,T') can be calculated in a similar way as in the Euclidean flat case. 
We have 

D,{r,r') = ^^6\k, ^0 (e-a'^^+z^^)!^-^'!) - e'^i^^^^^^^^'))^ , (72) 

or, in the "coordinate" space: 

D{t,t';x,x') = j dfl{k)uk{x)ul{x')Dk{T,T') = 

J dmn,(x)ul{x')^j^l-^^ (e-((^^+.^)l-^'l) - e-(('=^-^'^^)(^H-^'))) , (73) 

where / dfl{k) is given by Eq.(66). 

Now, let us apply this method for the case of a self-interacting theory with an interaction 
action given by: 

Sj[^] = jd'x^)j^<p\x). (74) 

In the same way, we can solve the equation using a perturbative series in A. The two-point function 
up to the one loop level is given by: 

( ^(n, Xi)(^(t2, X2) = (a) + {h) + (c), (75) 

where (a) is the zero order two-point function given by Eq.(72) and (b) and (c) are given respec- 
tively by 

iP) = -^5\krM) j dTGk,{ri-r)Dk{T,T)Dk,{T2,T), (76) 
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{c)^-^S\h,k2) J dmfj dTGk,{T2-T)Dk{T,T)Dk,{n,T), (77) 

where / d(3{q) = ^ / dkp J2q,j,M ^q^q- These are the contributions in first order. A simple 
computation shows that we obtain, at equal asymptotic Markov parameters (ti = r2 — > oo): 

Now we define the quantity / as: 

= /d*,i:X,-X,^. (79) 

where = A;^ + /U^. It is easy to show that the series in this equation is divergent. So, we need a 
procedure to regularize it and obtain a finite quantity for the two-point function. As we will see 
now, this can be done within the covariant stochastic regularization [36]. 
The generalization of Eq.(35) to general four-dimensional spacetimes is 

-^^{r,x) = --^ -^^Ux)=<p(r,x) + J d'^yy/gR^y(A)r](r,y), (80) 

where, for covariant reasons, the regulator is now a function of the covariant Laplacian: 

A,y = J (iS(V^),,(V^),„ (81) 

with 

iV,)xy = VlS\x-y). (82) 
Using the mode decomposition given by Eq.(59), the Eq.(80) reduces to 

d 1 

— <^fc(r) = --{e + + Vk{r)Rk, (83) 

where Rk = -R^j/(A)|A=-fe2+x and 

R^y{A) = I djl{k)uk{x)ul{y)Rk. (84) 
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Then, the solution for the regularized Langevin equation follows: 

ipkir) = r dr'Gkir, r')rj(r', k), (85) 
Jo 

where Gk{T,T') is given by Eq.(69). Then, it is easy to show that the zero order two-point 
function and the regularized contribution (6)a at the first order at the coupling constant are given 
by, respectively: 

D,{r, r') = j^^5\k, k')j^^^ (e-(('=^+M^)l-^'l) - e-((^^+-^)(^+^'))) , (86) 

and 

Let us isolate the part where we had problems. We have: 

= ^/*.E^X,-X,-^. (88) 

After some tedious calculations [47] [48] [49], we arrive at: 

= 1 e^+cx^^' dx{x^ - fi'^y/^erfciV^x), (89) 
4 J n 

where erfc{x) is the complementary error function. At the expense of considerable labor, it is 
possible to extend our results to the four-point function. 

A quite important point is that this regularization procedure preserves all the symmetries of 
the unregularized Lagrangian, since it is not an action regularization. The next step would be 
to isolate the parts that go to infinity in the limit A — > oo and remove them with a suitable 
redefinition of the constants of the theory, i.e., carry out the renormalization program. A natural 
question now would be if we can actually renormalize all the n-point functions at all orders at 
the coupling constant A. Birrel [50] has given arguments that a priori we cannot expect that a 
regularization independent proof of the renormalizability of the X(fi'^ theory in a curved background 
exists. One attempt of general proof of renormalizability of A<^^ theory defined in a spacetime 
which can be analytically continued to Euclidean situation was given by Bunch [51]. Using the 
Epstein-Glaser method, Brunetti and Fredenhagen [52] presented a perturbative construction of 
this theory on a smooth globally hyperbolic curved spacetime. 
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Our derivation shows that the stochastic regularization may be an attempt in a direction of 
such regularization independent proof, even though we are still restricted to the same situation 
studied by Bunch. Indeed, we know that, in this case of the Einstein universe, the free two-point 
function presents itself as a sum of Minkowski space positive frequency two-point functions, so it 
is natural to expect that, as it is in the Minkowski gularization independent proof of the 

renormalizability of such interacting theory does exist, at least for the Einstein universe. 

4 Stochastic quantization for the {X(p'^)d scalar theory: the 
Rindler case 

Now let us consider the stochastic quantization in a spacetime with an event horizon. The 
Schwarzschild black hole is the most familiar solution to the vacuum Einstein equations, which 
has an event horizon. At the origin, there is a curvature singularity and at i? = 2M we have 
a singularity due to the bad behavior of this particular coordinate system. Although it is not 
possible to extend the spacetime across the singularity, it is possible to analytically continue the 
manifold from r > 2M to the region r < 2M, (r 7^ 0). The maximal extension of the manifold 
described by the Schwarzschild line element with 2M < r < 00 is the usual Kruskal extension 
[53]. The Kruskal spacetime defines two outer asymptotically fiat regions and also two regions 
inside the event horizon, bounded by the past and future singularities (see Appendix A). Since 
we know that, close to the horizon, the Schwarzschild coordinates t and r behaves as Rindler's 
spacetime coordinates, we shall investigate the stochastic quantization in the Rindler metric. The 
background material relevant for us can be found in the text book [54] and also in the review 
articles [55] [56] [57] [58]. 

Let us consider a d-dimensional fiat Minkowski spacetime where we are using the usual cartesian 
coordinates — y^, y*^"^). It is possible to find a curvihnear coordinate system — 
{t,x^, ...,x'^~^,z), called the Rindler's coordinate system and can be shown that this coordinate 
system is one naturaly adapted to an observer with constant proper acceleration. It is important 
to stress that this coordinate system with the respective coordinate transformation cover only the 
region > \y^\. Since this coordinate system does not cover the whole Minkowski spacetime, 
we can define three coordinate transformations with the respective coordinate systems defined 
in different regions of the Minkowski spacetime. These regions are known in the literature as 
Rindler's L, Milne F, and finally Milne P. The four coordinate transformations and the coordinate 
systems together cover all the Minkowski spacetime. The coordinate system that cover the region 
inside the hght cone is the d-dimensional Milne spacetime. 

The implementation of the canonical quantization in Rindler's spacetime is very simple, since if 
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the spacetime has a stationary geometry, there is a time-hke kilhng vector field K that generates 
one-parametric Lie group of isometries and the orthogonal modes which satisfies CKUn{x) — 
—ioj Un, where Ck is the Lie derivative with respect to K and the Un{x) are the positive frequency 
modes. In this situation, there is a natural way to define positive and negative frequency modes. 
Note that the Rindler's line element is t independent, and consequently there is a straightforward 
way to define positive and negative frequency modes in order to impose the canonical quantization 
in Rindler's spacetime. 

It has long be recognized that the Rindler's vacuum | 0, R) that appear in the FuUing canonical 
quantization of a scalar field in the Rindler's spacetime is not unitarily equivalent to the Minkowski 
vacuum | 0, M >. Let us define the following generating functionals, Z{h) and ZR{h), i.e., 

^.^^ ^ (0,M, out\0,M, in)h .gQ. 

and 

_ (0,i?, out\Q,R,in)y^ 
^^^^>- {0,R,out\0,R,tn)^^, ^ 

where | 0, M, in), \ 0, M, out), \ 0, R, in) and finally | 0, R, out) are the IN and OUT vacuum states 
for the theory with the Minkowski and Rindler Hamiltonian density respectively. 

Starting from Z{h) it is possible after analytic continuation and imposing periodicity in the 
Euclidean time, to define Z{l3] h), i.e., the finite temperature Schwinger functional. It is clear that 
limj3^ooZ{f3; h) = Z{oo; h) = Z{h), where Z{oo; h) is the zero temperature Euclidean functional 
which generates the Schwinger functions in the whole Euclidean space. At the same way it is 
possible to define ZR{j3;h), i.e. the finite temperature Euclidean Rindler functional. Again we 
have limi3^ooZR{f3; h) — ZR{oo;h) = Z{h), where Zn{oo;h) is the zero temperature Euclidean 
functional which generates the Schwinger functions in the analytic extended Rindler manifold. 

From the above definitions, let us define the following Schwinger functions. We are following 
the discussion of Christensen and Duff [59]: 

i) Gp{x, x') two-point function, i.e. the Schwinger function obtained from the finite temperature 
Schwinger functional. 

ii) Goo{x,x') two-point function, i.e. the Schwinger function obtained from the zero tempera- 
ture Schwinger functional. 

Now, let us analyse the two-point function that we are able to construct in the Euclidean 
manifold. Remember that the Rindler Euclidean metric possesses a coordinate singularity at the 
origin where z — The point ^ = is a conical point of the Euclidean manifold but is a regular 
point if the Euclidean time is periodic. 

i) G'^^\x,x') two-point function which is periodic in the Euclidean Rindler's time. This is the 
usual scalar two-point Schwinger function. It is a function of the geodesic distance beween two 
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points in the manifold. Since all path around the origin are topologically equivalent, therefore 
G^^\x,x') should be periodic in Euclidean Rindler's time. Therefore these Green's functions 
can be generated by functional derivatives in the Z{P;h), i.e., the finite temperature Schwinger 
functional. 

ii) G'^'^\x,x') two-point function. This two-point function appears in the situation where the 
Euclidean space has a hole in the origin. In the case where the two points arc on the circular 
world-line it is possible to relate the two two-point functions G''^\x,x') and G^^\x,x'). This can 
be done using a simple trigonometric identity 

G('\ry)^ E G°(r,r' + ^^). (92) 

n=-oo ^ 

Since the finite temperature Schwinger function must be periodic in the Euclidean time with 
periodicity P, making the identification /5 = ^ it is clear that G^^\x,x') = Gp{x,x'). Also we 
have G'^^\x,x') = Goo{x,x'). Therefore these Green's functions can be generated by functional 
derivatives in the Z{oo;h) which is the zero temperature Euclidean functional which generates the 
Schwinger functions in the whole Euchdean space. 

Now, after a Wick rotation, we should apply the stochastic quantization for the four-dimensional 
Rindler space, which becomes a multiple connected manifold, with the Euclidean metric: 

ds'^ ^ z^dt^ + dxj + dxl + dz'^. (93) 

In the Euclidean Rindler space, the Langevin equation for the field (f{x) reads: 

^^{t,x) = -V9 j^Ux)=vir,x) + V{r,x), (94) 

where 5*0 is the Euclidean action for the scalar free fields. Notice that this expression is quite 
different from the Langevin equation for the Einstein metric, Eq.(43), since we have inverse power 
of the determinant of the metric multiplied by 5^^- So, with the Rindler metric given above, we 
have: 

d 

— V?(r, x) = -(-A + z'^m'^)ip{T, x) + ri{T, x), (95) 

where the operator A is defined by: 

A = d^ + z\dl+dl + dl) + zd,. (96) 

To proceed with the implementation of the stochastic quantization, we have to use the generalized 
Einstein relations, given by Eq.(46) and Eq.(47). 
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From the modes presented by references [60] [61], we can obtain the following Euclidean modes: 

1 / ■t- 
Uky{x,z) = —[psm\i{Tiv)\ e'^-^'Ki^ifiz), (97) 

where k.x — kot + kiXi + k2X2, = \Jk'^ -\- m? and Kn{x) is the Macdonald function. This are 
normalized such as: 

■n^ J z (i^ + i/')(smh(7ri/) smh(7rz/'))^/^ 

Now, we use the general mode decomposition, given by Eq.(52) and Eq.(53), with the following 
measure: 

j dfl{k) = J^" duj dk = du J dki J dk2j (99) 

Remembering Eq.(97) for the modes, we will have again that each mode coefficient obeys a 
Langevin equation of the form: 

d 

^</'A;i/(t) = -(/Cq + ^'^)Vkv{T) + r}ku{r). (100) 

The solution for Eq.(lOO), with the initial condition ^Pku{T)\T=o = 0, reads: 

vUr) = r dr'GUr, t'Ht', k), (101) 
Jo 

where 

GUt,t') = exp {-{kl + u'){t - r')) 9{t - r') (102) 

is the retarded Green function for the diffusion problem. Using the Einstein relations, we get that 
the mode coefficients for the random noise satisfies 

(%(t)), = (103) 

and 

( Vk{r)^k'iT') ), = 25{t - t')5{K k')5{v, u'), (104) 

where S{k, k') = 5{ko + A;[,)5(fc + k'). 

The two-point function D1{t,t') can be calculated in a similar way as in the Einstein case. 
We have 

Dlir, t') = 5{k, k')5{u, ^Ott^^ (e-(('=e+^^)l--'l) _ e-(('=a+'^^)(-+-'))) , (105) 

[Kq -j- ly ) ^ 
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or, in the "coordinate" space: 



2 roo r 

D^(x,x';t,t') = -— / diy dkvsmh.{'Kv)Kiy{ixz)Kiy{iJiz') 
4:71^ Jo J 

tie ( CC ) 

^_ L^(ik'o+'^')\r-r'\) _ . (106) 

kQ + ^ ^ 

In the hmit r = r' oo, and proceeding with similar steps as in reference [60], one may prove 
that we have obtained the usual result for the free two-point function, i.e., the Schwinger two-point 
function for a (straight-line) geodesic distance between the points x e x'. As is well known [59], 
this function is determined uniquely by the requirement that, in the absence of any holes in space, 
one identifies the "angle" t (which is the Euclidean "time" ) with the angle t + 27m {n ^integer) 
because all paths around the origin are topologically equivalent. Thus, as expected, should be 
periodic in t with period 27r. But we know that, when we do the analytical continuation t ^ it in 
the Rindler metric, the event horizon, which in Minkwoski space is represented as an impenetrable 
barrier to the accelerating observer, translates in Euclidean language into the statement that closed 
paths around the origin cannot be continuously shrunk to a point. An inertial observer sees a 
space with the usual topology, R"^, but the topology seen by the accelerating observer is different, 
i.e., his space has a hole in it. Then, paths winding around the origin cannot all be deformed into 
each other; they fall into topologically distinct classes labelled by the winding number n. Now, 
we do not identify t with t + 27rn. This two-point function, which we shall call D^, exhibits an 
infinite periodicity. To obtain it, we should work with the following modes [62]: 

u^x, z) = 2^(i/sinh(7ri/))^/'e'(*^i^^+^^^^)-'=°l*lKi,(//^), (107) 
With this prescription, we obtain the same results as found in literature [63]: 

D%x,x';t,t') = D\x,x';t,t') 

n=—oo 

I n=+oo oo r 

— V / du dku&mh.{Tiv)KiJiiz)KiJiiz') 



n=—(X) 
J,{k\xi-\-k2X2)—ko\t-\-2'Kn\ 



^ (e-«'=o+'^')l---'l) - e-«^o+'^')(-+-'))) , (108) 



where the geodesic distance between the points x and x' is an arc length in the function D^. It is 
interesting to note that, using the modes given by Eq.(107) (see reference [62]), and then summing 
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up all configurations with a winding number n we reach the same result as using the modes given 
by Eq.(97) (see reference [60]). 

Let us repeat the method for the case of self-interacting theory with the interaction action 
given by: 

Jd'xy/^)^^\x). (109) 
In this case, we have similar equations to the Einstein case. Eq.(32) reads, in the Rindler space: 

(b)^~6{h,h)S{iy,,iy2) J dm £ dTG,,,,(n-T)Dl{T,T)Dl,^{T2,T), (110) 
which, in the limit r = r' — > oo, becomes: 

= "^^ (fci!.i)(fef+fci!.,^+.i) / ^^^^^ (FT^ ^'''^ 

with similar divergences and J d[3 = J du J dku smh.{'n-u)Ki,y{fiz)Ki,,{fiz). Now, we may apply 
the continuum regularization as before [36], with a regulator that is a function of the following 
operator: 

A' = z'{^^^ + dl + dl) + zd^ - zV. (112) 
The regularized Langevin equation reads: 

-Q:^fkv{r) = -{kl + v'^)ipkv{r) + RkuVkuir), (113) 

with the solution: ^ 

Vku{r)= f dT'Gkuir,T')RkuVir',k). (114) 

JO 

where Rkv{^) = i?(A'; A)| A'=iy2. Then, it is easy to show that the zero order two-point function 
and the regularized contribution (6)a at the one loop level are given by, respectively: 

DUr, r') = 6ik, k')6{., u')-^^ (e-((^^+'^^)l-'l) - e-«'=o^+^^)(^+^'))) , (115) 

and 

Wl~ = -^fe. few-.. (,f hi) w + . I + "'''^"'^wh- ("'^) 
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So, the regularized part of the contribution above reads: 



^ ^ An^ Jo Jo J '^^^^^H'^^)^i''if^^)^i''il^^) P + ^'^'^^^ 
With some tedious manipulations [47] [48] [49], we arrive at: 

where the function /(A, z) is given by 

„ f+°° da ( iisAi / A \ .HTsAi / A \ 

J-00 * y ^ er/c(^-^-^(zQ;+7r)J -e 4 er/c(^-^^(zQ;-7r)J 

where 7^ = 2z'^{l + cosh a) and erfc{x) is the complementary error function, which satisfies the 
identity: 

erfc{-x) = 2 - erfc{x). (120) 

It is not possible to present the solution of the integral in Eq.(119) in terms of known functions. 
In spite of this inconvenience, it is easy to see that this function /(A, z) has strong convergence, 
proving that we indeed regularize the contribution for the two-point function at the one-loop level. 
Similar calculations can be carried out for the four-point function. Prom the discussions it should 
be clear that is a simple matter to derive the Schwinger functions of the theory at least in the 
one-loop level. 

There are two rather subtle points that we have to investigate. The first one is related to the 
behavior of the random noise near event horizon, and the second one is the validity of the program 
beyond the Euclidean signature. 



5 Conclusions and perspectives 

The stochastic quantization method was used to study self-interacting fields in manifolds which 
can be analytically continued to the Euclidean situation, i.e., static Ricmannian manifolds, namely, 
the Einstein manifold and the Rindlcr manifold. First, we have solved a Langcvin equation for 
the mode coefficients of the field, then we exhibit the two-point function at the one-loop level. It 
was shown that it diverges and we have used a covariant stochastic regularization to regularize it. 
It was shown that, indeed, the two-point function is regularized. 
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A natural question that arises, when we work in Rindler space, is that what happens to the 
noise field correlation function, given by Eq.(47), near an horizon. From this equation, we see that, 
whenever we have g — det^^j, = 0, this correlation function diverges, and, therefore, all n-point 
correlation functions ( v?(ri, xi)ip{T2, X2).-.'^{Tn, Xn) )r, will have meaningless values, in virtue of the 
solution of the Langevin equation. We may implement a brick wall-like model [64] [65] in order to 
account for these effects; in other words, we may impose a boundary condition on solutions of the 
Langevin equation at a point near the horizon. On the other hand, in the limit g = det g^i, oo, 
all the correlation functions vanish. We would like to point out that another and powerful way 
to understand the behavior of the correlation functions in the limit g — det g^^ = may be the 
stochastic quantization of the two-dimensional Schwinger model in a curved background. This 
can also be a way to study the most pertinent form of the Langevin equation to be used in curved 
spacetimes. This subject and also the link between our results and the program to implement the 
perturbative renormalization in domains where the translational invariance is broken [66] [67] [68] 
[69] is under investigation by the authors. 

We discussed also that the presence of the Markov parameter as an extra dimension lead us 
to a regularization scheme, which preserves all the symmetries of the theory under study. This is 
a quite important point in a regularization procedure. It must preserve all the symmetries of the 
unregularized Lagrangian. Although many authors have stressed that we can not expect that a 
regularization independent proof of the renormalization of theories in a curved background exists, 
since the stochastic regularization is not an action regularization, may be a way to construct such 
proof. We are aware of the fact that the stochastic quantization program can be implemented 
without problems (modulo event horizon problems, etc), if it is possible to perform the Wick 
rotation, obtaining a real Euclidean action. 

The picture that emerges from the discussions of the paper is that the implementation of the 
stochastic quantization in curved background is related to the following fact. For static manifold, it 
is possible to perform a Wick rotation, i.e., analytically extend the pseudo-Ricmannian manifold to 
the Riemannian domain without problem. Nevertheless, for non-static curved manifolds we have to 
extend the formalism beyond the Euclidean signature, i.e., to formulate the stochastic quantization 
in pseudo-Riemannian manifold, not in the Riemannian space (as in the original Euclidean space) 
as was originally formulated. See for example the Refs. [11] [12]. In the Ref. [11], the authors 
proposed a modification of the original Parisi-Wu scheme, introducing a complex drift term in the 
Langevin equation, to implement the stochastic quantization in Minkowski spacetime. Gozzi [12] 
studied the spectrum of the non-self-adjoint Fokker-Planck Hamiltian to justify this program. See 
also the Refs. [70] [71]. Of course, this situation is a special case of ordinary Euclidean formulation 
for systems with complex actions. 

The main difference between the original Euclidean and the implementation of the stochastic 
quantization in Minkowski spacetime is the fact that in the usual case a real Euclidean action the 
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approach to the equiUbrium state is a stationary solution of the Focker-Planck equation. In the 
Minkowski formulation, the Hamiltonian is non-Hermitian and the eigenvalues of such Hamiltonian 
are in general complex value. The real part of such eigenvalues are important to the asymptotic 
behavior at large Markov time. The approach to the equilibrium is achieved only if we can 
show its positive semi-definiteness. The fundamental question is: what happens if the Langevin 
equation describes diffusion around complex action? Some authors claim that it is possible to 
obtain meaningful results out of Langevin equation diffusion processes around complex action. 
See for example the seminal paper of Parisi [72]. Klauder et al. [73] investigated the complex 
Lagevin equation, where some numerical simulations in one-dimensional systems was presented. 
See also the Refs. [74] [75] . Finally we would like to mention the approach developed by Okamoto 
et al. [76], where the role of the kernel in the complex Langevin equation was studied. 

We would like to remark that there are many examples where Euclidean action is complex. The 
stochastic quantization in Minkowski spacetime, as we discussed; system with chemical potential 
as for example QCD with non-vanishing chemical potential at finite temperature; for SU{N) 
theories with N > 2, the fermion determinant becomes complex and also the effective action. 
Complex terms can also appears in the Langevin equation for fermions, but a suitable kernel 
can circumvent this problem. Another quite instructive case that deserves our attention is the 
stochastic quantization of topological field theory. One of the peculiar feature within this kind 
of theory is the appearance of a factor i in front of the topological action in Euclidean space. 
In a topological theory, the path integral measure weighing remains to be e*'^, even after the 
Wick rotation. An attempt to use a Markovian Langevin equation with a white noise to quantize 
the theory, fails since the Langevin equation will not tend to any equilibrium at large Markov 
parameter. In the literature there are different proposed to solve the above mentioned problem. 
In a pure topological Chern-Simons theory, Ferrari et al. [77] introduced a non-trivial kernel 
in the Langevin equation. Other approach was developed by the Menezes et al. [28]. These 
authors showed that using a non-Markovian Langevin equation with a colored random noise, the 
convergence problem can be solved. These authors proved that it is possible to obtain convergence 
toward equilibrium even with an imaginary Chern-Simons coefficient. We conclude saying that 
several alternative methods have been proposes to deal with interesting physical systems where the 
Euclidean actio is complex. These methods do not suggest any general way of solve the particular 
difficulties that arise in each situation. 

The program of using the complex Langevin equation for an investigation of quantum field 
theories in Minkowski spacetime is not yet established and is still under discussion. Its clear that 
until the problem that we discussed above is solved, one can not say for sure that the imple- 
mentation of the program of stochastic stochastic quantization in non-static pseudo-Riemannian 
manifolds can be successfully implemented. 
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A The Kruskal extension of the Schwarzschild spacetime. 

The spacetime of a spherically symmetric non-rotating and eletrically neutral body is described 
by the Schwarzschild metric 



where dfl^ = dO^ -\- siv? 6 d(f) is the metric of the unit two-sphere. The coordinate system (t, r, 6*, (f)) 
provides a frame in which the metric components are time-independent. Since, it is possible to 
find a Killing vector field which is hypcrsurfacc orthogonal to the family of spacelikc hypcrsurfaces 
t = cte, the Schwarzschild solution is also static. The line element given by Eq.(A.l) describes the 
external gravitational field generated by any spherical mass, whetever its radius [78] . This includes 
several interesting situations, as a spherically symmetric star which undergoes radial pulsations 
or ever a radial spherically symmetric gravitational collapse. In the extreme case of a complete 
gravitational collapse, we have to consider the Schwarzschild line element as an empty spacetime 
solution for all values of r. Since the line element has two singularities, one at r = and another 
at r = 2M, it represents only one of the patches < r < 2M or 2M < r < oo. It is not difficult 
to show that although the Schwarzschild line element is singular at r = 2M, all the invariants 
constructed with the Riemann tensor and its contractions are well-behaved at r = 2M, and this 
point is a singularity due to an inappropriate choice of coordinates. 

On the other hand, choosing the Schwarzschild line element to describes the patch < r < 2M, 
we have that the curvature scalar diverges at r = 0, since we have R^^p^R^'^P'' = Therefore 
this point is a real spacetime singularity. Thus, although it is not possible to extend the spacetime 
across the singularity, it is posssible to analytically continue the manifold from r > 2M to the 
region r < 2M, (r 7^ 0). The maximal extension of the manifold described by the Schwarzschild 
line element with 2M < r < 00 was obtained by Kruskal. 

From the condition for null geodesies it is possible to show that the radial null geodesies in the 
Schwarzschild spacetime are given by 



(CNPq). 




(A.l) 




(A.2) 
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and allow us to define r* as 



* =r + 2Mln(^-l), (A.3) 



and it is clear that the radial geodesies must satisfy t — r* + cte or t — —r* + cte. 

The next step is to define the null coordinates u and v, where u — t — r* and v — t -\- r*. 
Obviously, it is possible to write the Schwarzschild line element in terms of the null coordinates u 
and V, and we have 

c/s^ = e 4M dudv + r'^dQ.'^, (A.4) 

r 

where the metric of the unit two-sphere is inalterated. To go further, let us define new coordinates 
U and V, where 

f = exp(^), (A.5) 



and 



V^ = exp(— ), (A.6) 



Using the coordinates U and V , the metric can be written as 

ds^^- ^^'^"'~^ dUdV + rHn\ (A.7) 
r 

The value r = 2M is no more a singularity, since corresponds to U = or K = 0. Finally, to 
obtain the Schwarzschild metric in the Kruskal form we have only to define the coordinates T and 
X, by choosing T = ^{U + V) and X = — U). The final form of the Schwarzschild metric in 
terms of the Kruskal coordinates (T, X.O^cf)) is 

ds-" = 32M3e-2M ^ ^^2) ^ ^2^j^2_ (A.8) 

r 

Note that the coordinate transformation between the original coordinates {t, r) and the Kruskal 
coordinates (T, X) is given by 

^'-^'=f;^7-l')expf-^V (A.9) 



V2M / ^\2M 

and 

t = 4Mtanh-^ (^^) . (A.IO) 
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Therefore, we show how it is possible to analytically continue the manifold from r > 2M 
to the region r < 2M, and the event horizon, a sphere of radius r = 2M is only a coordinate 
singularity, which can be removed by a suitable coordinate transformation. Althought the apparent 
singularity at the horizon has disappeared, there are true singularities in the Kruskal extension 
of the Schwarzschild spacetime. The physical singularity at r = corresponds to the values X = 
\/T'^ — \ and X = —\/T'^ — 1. The original Schwarzschild solution for r > 2M corresponds to the 
region where observers can obtain information. Since there are two event horizon, the future event 
horizon and the past event horizon, the Kruskal spacetime defines two outer asymptotically flat 
regions and also two regions inside the event horizon, bounded by the past and future singularities. 
The black-hole is the maximal analytic extension of the vacuum Schwarzschild solution. It is well 
known that close to the horizon, the Schwarzschild coordinates t and r behaves as Rindler's 
spacetime coordinates. 
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